Spin Dynamics in the Magnetoelectric Effect LiCoP04 Compound 
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Inelastic neutron scattering (INS) experiments were performed to investigate the spin dynamics in 
magnetoelectric effect (ME) LiCoPCU single crystals. Weak dispersion was detected in the magnetic 
excitation spectra along the three principal crystallographic axes measured around the (0 1 0) 
magnetic reflection. Analysis of the data using linear spin-wave theory indicate that single-ion 
anisotropy in LiCoPCU is as important as the strongest nearest-neighbor exchange coupling. Our 
results suggest that Co 2+ single-ion anisotropy plays an important role in the spin dynamics of 
LiCoPCU and must be taken into account in understanding its physical properties. High resolution 
INS measurements reveal an anomalous low energy excitation that we hypothesize may be related 
to the magnetoelectric effect of LiCoPC>4. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.30.Et 



I. INTRODUCTION 



LiCoP04 is an antiferromagnetic (AFM) insulator be- 
longing to the olivine family of lithium orthophosphates 
that share the general chemical formula LiMPC>4 (M 
= Mn 2+ , Fe 2+ , Co 2+ , Ni 2+ ) with four formula units 
per unit celli«2. These materials continue to attract 
much attention due to their exceptionally large magne- 
toelectric (ME) effect and the anomalies exhibited in 
the ME coefficients as a function of temperature and 
magnetic field 3 ' 4 '-. To date, it remained an open ques- 
tion whether the anomalies observed in the ME effect of 
LiMP04 are intrinsic due to the particular local environ- 
ment surrounding the transition metal ions in LiMP04 
or due to domain formation structure. The local envi- 
ronment can be slightly distorted in a magnetic (elec- 
tric) field by virtue of the spin-orbit coupling, giving 
rise to a collective ferroelectric response (magneto-ferro- 
elastic effect The recent domain structure observed 
by second harmonic generation (SHG) in LiCoP04 was 
attributed to coexisting AFM and ferrotoroidic domains 
which may play a role in giving rise to the ME effect 8 . 
In addition to the strong ME effect, this family of ma- 
terials also exhibits intriguing magnetic properties. At 
low temperatures, L1MPO4 systems undergo transitions 
to AFM long-range-order (LRO), adopting similar mag- 
netic structures, differing only in spin orientation. For 
example, with increasing temperature, LiNiPC^ first un- 
dergoes a first-order commensurate-incommensurate (C- 
IC) phase transition at Tm w 20.8 K changing from a 
co-linear AFM state to a long-range IC order state, fol- 
lowed by a second-order phase transition from long-range 
IC to short-range IC order at T IC w 21.7 K^. On the 
other hand, LiMnPC>4 undergoes an AFM LRO tran- 
sition at Tn ~ 34 K as well as a field induced spin- 
flop transition^. Weak ferromagnetismii, an ME "but- 
terfly loop" anomaly^, and strong magnetic anisotropy 
have been observed in LiCoP04 13 which also exhibits 
the largest ME coefficient 5 among its counterpart com- 
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FIG. 1: (Color online) Magnetic unit cell of LiCoP04 dis- 
placed (0.25 0.25 0) r.l.u compared to the atomic unit cell. 
Only the Co 2+ magnetic ions are shown for clarity. The intra- 
plane (solid) and inter-plane (dashed) magnetic exchange in- 
teractions considered in the spin Hamiltonian (Eq. [1} are la- 
beled. It is assumed that the spin is oriented strictly along 
the crystallographic 6-axis in the model calculation. 



pounds, making it of particular interest to study. 

LiCoP04 crystallizes in an orthorhombic symmetry, 
space group Prima (no. 62) at room temperature, with 
lattice parameters a = 10.093, b — 5.89, and c = 
4.705 A. The structure consists of buckled CoO layers 
stacked along the crystallographic a-axis and the mag- 
netic Co 2+ (S = 3/2) ions are surrounded by oxygen 
ions in a strongly distorted C0O6 octahedral coordina- 
tion. LiCoP0 4 develops AFM LRO at T N w 21.8 Ki 2 .^. 
Earlier studies have indicated a simple two-sublatticc 
AFM state below Tjv with spins aligned along the 6-axis 
(whereas LiFeP04, LiMnP04, and LiNiP04 have spins 
oriented along the b, a, and c-axis respectively ^il 5 .. Fig. 
Q] illustrates the magnetic structure of LiCoP04. For 
simplicity, only the Co 2+ magnetic ions are shown. As 
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depicted in Fig. [T] different magnetic exchange pathways 
are at play in this compound. Within the buckled CoO 
layer, nearest-neighbor Co 2+ ions are strongly coupled 
( Ji) through Co-O-Co super-exchange interactions, while 
additional inplane magnetic interactions between next- 
nearest- neighbors ( J2, J3) are mediated through the PO4 
group. Between adjacent layers, the magnetic couplings 
of nearest-neighbor ions ( J4, J5) are also mediated via the 
PO4 groups. The PO4 coupling are found to be rather 
strong and cannot be neglected^^. Due to its layered 
structure, LiCoPC>4 exhibits properties between two- and 
three-dimensional (2D and 3D) magnetic systems. 

Although the physical properties of LiCoPCU have 
been studied extensively in the past, there remain 
a number of puzzles in this compound^. Re- 
cent magnetoelectric&i2, magneto-optic^, and magnetic 
property^ studies do not agree with the originally pro- 
posed co-linear AFM structure and suggest a more com- 
plex magnetic structure for LiCoPC>4. Neutron diffrac- 
tion studies suggest the moments in LiCoPC>4 in the 
AFM phase are not strictly aligned along the 6-axis 
but are uniformly rotated from this axis by a small 
angle (~ 4.6°)—. The observation of weak ferromag- 
netism and an ME "butterfly loop" anomaly further 
motivated studies of LiCoPC>4 both experimentally and 
theoreticall y 20 ' 21 ^ 22 . In this paper, we report single crys- 
tal inelastic neutron scattering (INS) studies that yield 
the microscopic magnetic interactions in LiCoPCU. The 
data were analyzed within the linear spin wave approx- 
imation using a spin Hamiltonian explicitly including 
the intra-, inter-plane nearest neighbor, next-nearest- 
neighbor exchange interactions and single-ion anisotropy, 
which are determined in this study. 



II. EXPERIMENTAL TECHNIQUES 

All measurements reported here were carried out on 
LiCoPC>4 single crystals. Large crystals were grown for 
INS experiments by a LiCl flux method similar to that 
reported in Ref. |23| . High purity starting materials 
of CoCl 2 (99.999%, MV laboratory), Li 3 P0 4 (99.999% 
Aldrich) and LiCl (99.999%) were thoroughly ground to- 
gether at a molar ratio of Li3P04:CoCl2:LiCl = 1:1:1 and 
sealed in a Pt crucible under Ar atmosphere. Note that 
the reaction of Li 3 P0 4 + CoCl 2 + LiCl — ► LiCoP0 4 + 
3LiCl yields a molar ratio of 1:3 between LiCoPC>4 and 
the flux material LiCl which we found crucial in grow- 
ing large LiCoP04 single crystals. The mixed powder is 
pre-melted at 800° C and slowly heated to 900° C. The 
crucible was maintained at 900°C for 10 hours and then 
slowly cooled to 640°C at a cooling rate of 0.7°C/hr 
and then furnace-cooled to room temperature. Large 
purple color LiCoP04 single crystals were obtained and 
extracted by dissolving LiCl in water. The crystals 
were characterized by X-ray diffraction measurements 
and found to be of pure single phase. The lattice pa- 
rameters determined from our neutron diffraction mea- 



surements a = 10.159, b = 5.9, and c = 4.70 A at 8 K 
are in good agreement with prior results 14 . 

Two LiCoPC>4 single crystals grown from the same 
batch were used for the INS experiments. Sample #1, 
m w 0.8 g, and sample #2, m ss 0.4 g, were oriented in 
the (H K 0) and (0 K L) scattering plane, respectively. 
INS experiments were performed using the HB1A triple- 
axis spectrometer (TAS) at the High Flux Isotope Reac- 
tor (HFIR) neutron scattering facility at Oak Ridge Na- 
tional Laboratory, the BT7 thermal TAS, and the SPINS 
cold TAS at the NIST Center for Neutron Research at 
the National Institute of Standards Technology (NCNR) . 
The magnetic excitations along the (H 1 0) and (0 K 0) 
directions were measured using the HB1A TAS on sample 
#1 which was mounted on a thin aluminum disk, sealed 
in an aluminum sample can and kept under helium at- 
mosphere and cooled using a closed-cycle He refrigera- 
tor. Collimations of 48'-48'-sampZe-40'-68' downstream 
from reactor to detector were used throughout the ex- 
periment. Constant wave- vector scans were performed 
at T = 8 K (fully ordered state) and T = 35 K (well 
above Tjy). Spin excitations along the (0 1 L) direction 
were measured using the BT7 TAS at NIST on sample 
#2. A fixed final energy of Ef = 14.7 meV and a open- 
bO'-sample-SO'-open collimation were used with pyrolytic 
graphite (PG(0 2)) analyzer crystals in flat mode. High 
resolution INS measurements were carried out on sample 
#1 using the SPINS TAS with a fixed final energy, E f 
= 5 meV. A collimation of open-80-sample-80'-open was 
used with a cold Be filter in the scattering beam. Polar- 
ized neutron scattering experiments were also performed 
using the BT7 TAS with a fixed final energy of Ef = 
13.7 meV to clarify the nature of the ~ 1.2 meV low en- 
ergy excitation. The polarization analysis technique as 
applied to this study is discussed in Ref. I2ll25fl26l 3 He 
spin filters (polarizer) are mounted before and after the 
sample with a spin flipper in the incident beam. The 
sample is maintained in a horizontal or vertical magnetic 
guide field of ^5 Oe such that the neutron polarization 
p is parallel to the momentum transfer Q, p || Q when a 
horizontal field is applied at the sample position, or p_LQ 
when a vertical magnetic field is applied. With the spin 
flipper off, we measure the (-1 — h) non-spin-flip scattering. 
On the other hand, with the spin flipper on, the (-+) 
spin-flip scattering is measured. Sample #2 oriented in 
the (0 K L) scattering plane and a open-HO'-sampleSO 1 - 
open collimation were used throughout the polarization 
measurements. Constant wave-vector scans were carried 
out with both (++) and (-+) configurations. All mea- 
surement results have been normalized to a beam monitor 
count. 



III. EXPERIMENTAL RESULTS AND 
MODELING 

Representative constant wave- vector scans with energy 
transfer between Klu = 2.5 meV and 8 meV measured us- 
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FIG. 2: (Color online) Representative constant wave-vector 
scans plotted as scattering intensity versus energy transfer 
between hui = 2.5 meV and 8 meV. (a). Temperature depen- 
dence of the ~ 4.7 meV excitation measured at (0 1 0) at T 
= 8 K and 35 K using the HB1A TAS. (b). T = 8 K constant 
wave-vector scans using the BT7 TAS measured at (0 1 0) 
and (0 1 1.5) showing a single magnetic excitation that ex- 
hibits modest but unambiguous dispersion along the (0 1 L) 
direction. Intensities were normalized to the incident neutron 
flux by counting against neutron monitor counts. 



ing the HB1A and BT7 TAS are shown in Fig. H Note 
that error bars in this paper are statistical in origin and 
represent one standard deviation. Fig. [2] (a) shows the 
temperature dependence of the magnetic excitation mea- 
sured at (0 1 0). At T = 8 K in the fully ordered phase, 
a single excitation with energy transfer of hu> ~ 4.7 meV 
is detected. At a temperature well above TV (T = 35 K) 
the peak disappears demonstrating the excitation is mag- 
netic in origin. Below TV, further detailed measurements 
at (0 1 0) as a function of temperature indicate no signif- 
icant temperature dependence of this excitation. Fig. [2] 
(b) shows the T = 8 K constant wave-vector scans mea- 
sured at (0 1 0) and (0 1 1.5), which typically correspond 
to the minimum and maximum spin wave excitations. 
The data clearly show that the excitation observed at 
hu) w 4.7 meV at (0 1 0) shifts to higher energy transfer 
Tiuj ~ 5.3 meV at (0 1 1.5), yielding a maximum energy 
shift of ~ 0.6 meV. This indicates that the overall disper- 



sion along the (OIL) direction is modest compared to 
an exchange energy of fcTjy ~ 2 meV. Similar behaviors 
were observed for the dispersions along the (0 K 0) and 
(H 1 0) directions. Note that data in Fig. [5] (a) shows 
strong scattering at 2.5 meV, the lowest data point plot- 
ted. We will show later that it is from scattering of a low 
energy excitation at ~ 1.2 meV which has been observed 
in recent SPINS high resolution measurements. 

To determine the magnetic excitation spectra along all 
three principal axes directions, a series of constant wave- 
vector scans were carried out in the (H K 0) and (0 K 
L) scattering planes at T — 8 K around the (0 1 0) mag- 
netic reflection. Fig. [3]depicts the ground state magnetic 
dispersion relations along the (H 1 0), (0 K 0), and (0 1 
L) directions constructed from energy scans at constant 
wave- vector. We determined the peak positions assuming 
Gaussian peak-shapes that were fit to each of the con- 
stant wave-vector scans measured. For both the HB1A 
and BT7 triplc-axis-spectrometers, the energy resolution 
at the elastic position was AE w 1 meV. The experimen- 
tal uncertainties had a significant effect in the theoretical 
modeling as described in the text below. The measured 
spectra indicate a spin wave excitation of hui ~ 4.7 ±0.24 
meV at (0 1 0) which vanishes abruptly above T/v, while 
modest dispersion was observed along all three principal 
symmetry directions, with the scale in Fig. [3] being cho- 
sen to best exhibit the dispersion that falls within a band 
of 0.8 meV. This relatively weak dispersion suggests an 
Ising-like model in LiCoPCU; in a pure Ising model the 
magnetic excitations are completely dispersionless. 

To analyze the measured spin wave dispersion curves 
of LiCoP04 using linear spin wave theory, we consider 
the different magnetic exchange interactions as illus- 
trated in Fig. Q] and assume an AFM ground state with 
spins pointing strictly along the &-axis. Taking into ac- 
count the intra-plane and inter-plane nearest-neighbor, 
next-nearest-neighbor interactions, and the single-ion 
anisotropy, the Spin Hamiltonian can be expressed in the 
following form: 

i.j i,a 

where D a (a = x,y,z) represents the single-ion 
anisotropy along the x , y, and z-directions. In order 
to have the spins pointing along the z-axis in the model 
calculation, the Cartesian coordinate x, y, and z direc- 
tions are defined to be along the crystallographic a, c, 
and b directions respectively. The zero point of the en- 
ergy spectrum is chosen such that D Z =Q. Within a linear 
spin wave approximation, the derived spin wave disper- 
sion from Eq. |T| is given by: 

hu = yjA 2 - (B±C) 2 . (2) 

where 

A = 4S*(J! + J 5 ) - 25[J 2 (1 - cos(q • r 5 )) + J 3 (l - cos(q • 
r 6 )) + J 4 (2 - cos(q ■ r 7 ) - cos(q • r 8 ))] + D X S + D y S, 
B = D X S - D y S, 
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FIG. 3: (Color online) Spin-wave dispersion curves along three reciprocal directions constructed from a series constant wave- 
vector scans measured at T — 8 K. Data points are obtained based upon a Gaussian peak approximation. The solid and dash 
lines are calculations based upon a global fit to the linear spin wave approximation theory as described in the text. 



C = 2J 1 S(cos(c[ • ri) + cos(q ■ r 2 )) + 2J 5 S(cos(q ■ r 3 ) + 
cos(q • r 4 )). 

and r 4 denotes the vectors directed between two Co 2+ 
ions 

r x = (0, ft/2, c/2),ra = (0, 6/2, -c/2), 
r 3 = (a/2,6/2,0), r 4 = (a/2, -6/2,0), 
r 5 = (0,6,0), r 6 = (0,0, c), 
r 7 = (a/2, 0, c/2), r 8 = (a/2, 0, -c/2). 

Non-linear-least-squares fits of the spin-wave disper- 
sion expressed by Eq. [2] to the observed magnetic spec- 
tra yields: Ji = 0.743 ± 0.187 meV, J 2 = 0.105± 0.159 
meV, J 3 = 0.194± 0.131 meV, J 4 = -0.163± 0.08 meV, 
J 5 = -0.181± 0.125 meV, D x = 0.718± 0.192 meV, 
D y = 0.802± 0.208 meV. The obtained microscopic in- 
teraction parameter Ji is significantly larger than J 4 and 
J5 consistent with previous observations that the mag- 
netic behavior of LiCoP0 4 is intermediate between a 2D 
and 3D system. Moreover, the obtained positive value 
of Ji indicates in-plane AFM nearest-neighbor coupling, 
whereas negative J 4 and J5 values suggest inter-planc 
FM coupling along the a-axis consistent with the mag- 
netic structure, where the magnetic unit cell is doubled 
along the b- and c-axes, but not along the a-axis. J 2 
and J3 have the same sign as J\ indicating that they 
compete with J\ and may cause frustration, however, 
they are relatively weak compared to J\ (J2/J1 ~ 0.14, 
J3/J1 rs 0.26). Both D x and D y are positive favoring a 
ground state with the magnetic moment along the 6-axis 
consistent with the elastic magnetic neutron scattering 
results. An important result in our study is the large val- 
ues of the single-ion anisotropy compared to the nearest- 
neighbor coupling D x ~ D y ~ J x , Although strong 
single-ion anisotropy in LiCoP0 4 has been suggested by 
several models&^l, this study provides experimental ev- 
idence that the single-ion anisotropy is as important as 
the strongest magnetic exchange interaction in LiCoP0 4 . 
Such relatively strong anisotropy may split the S = 3/2 
quartet of the Co 2+ ion into two doublets rendering the 
suggested Ising-like character to LiCoPO^i^. 

The "±" sign in Eq. © (B ± C) indicates that there 



are two non-degenerate spin wave branches which come 
directly from the different values of D x and D y . Using 
the obtained best-fit parameters, the calculated disper- 
sion curves of the two branches are plotted as solid lines 
("B - C" branch) and dash lines ("B + C" branch ) in 
Fig. [3] The two calculated spin wave branches predict a 
maximum separation of ~ 0.3 meV at (0 1 0), (1 1 0), 
and (0 1 1). In the thermal neutron TAS measurements 
using HB1A and BT7 with a resolution of ~ 1 meV, 
only one excitation was observed at these wave vectors. 
We have two high resolution measurements using SPINS 
TAS, fixed Ef = 5 meV with a resolution of ~ 0.28 meV, 
at (1 1 0) and (0 1 0) with energy transfer up to 8 meV. 
The constant wave vector scan at (0 1 0) at T — 9 K is 
shown in Fig. 2] (a). At (0 10), where the model predicts 
the maximum separation between these two branches, 
only one excitation around ~ 4.7 meV is observed. The 
additional low energy excitation observed at ~ 1.2 meV 
does not agree with the model and is discussed below. 
Our results could not resolve the two branches for two 
possible reasons. First, the second excitation may be 
very weak in intensity, and our model does not predict in- 
tensities. Second, the intrinsic linewidth of the observed 
excitations are broader than the resolution (~ 1 meV) 
suggestive of contributions from both branches overlap 
and cause the broadening. 

High resolution measurements on SPINS show an 
anomalous low energy excitation below T/v that does not 
fit in the linear spin wave model. Fig. [4] (a) depicts the 
T = 9 K, 15 K, and 21 K data measured at (0 10). 
In addition to the Tito w 4.7 meV excitation, the T = 
9 K SPINS data clearly show a low energy excitation 
at huj w 1.2 meV. The peak position of this excitation 
is practically temperature independent but the peak in- 
tensity decreases with increasing temperature and even- 
tually the peak vanishes above Tjy. It also shows weak 
dispersion along all three reciprocal directions. Fig.[4](b) 
shows constant wave- vector scans measured along the (0 
K 0) direction at 9 K. Very weak dispersion was observed 
along this direction and the data show that this excita- 



CO 



C 
=J 

o 
O 



150 



100 



T = 9 K 
(0 1 0) 
(0 1.1 0) 
(0 1.2 0) 
(0 1.3 0) 




FIG. 4: (Color online) SPINS high resolution measurements 
of LiCoPCU. (a) Constant wave- vector scans measured at (0 1 
0)atT = 9K, 15K, and 21 K indicating a second low energy 
excitation at Tiu> ~ 1.2 meV below Tat. (b) Weak dispersion 
observed in the constant wave-vector scans measured along 
the (0 K 0) direction at T = 9 K. 



tion weakens in intensity (significantly) with increasing 
K and could not be detected at large K. Similar results 
were obtained along the (H 1 0) and (0 1 L) directions. 
In order to clarify the origin of this excitation, polarized 
neutron scattering experiments were carried out using 
the BT7 TAS. As discussed in Ref. HHHIH, coherent 
nuclear scattering is always non-spin-flip scattering (H — (-) 
because it never causes a reversal, or spin flip, of the neu- 
tron spin direction upon scattering. On the other hand, 
magnetic scattering depends on the relative orientation 
of the neutron polarization p and the scattering vector 
Q. Only those spin components which are perpendic- 
ular to the scattering vector are effective. Thus for a 
fully polarized neutron beam, with the horizontal field 
configuration, p || Q, all magnetic scattering is spin- flip 
scattering (-+), and ideally no non-spin-flip scattering 
will be observed. Our polarized measurements were car- 
ried out by performing constant wave-vector scans at (0 
1.2 0) at T = 7 K with (++) and (-+) configurations 
and a horizontal magnetic guide field at the sample po- 
sition (p || Q). Inset in Fig. [5] (a) first compares the (0 
2 0) nuclear Bragg scattering measured in non-spin-flip 
(H — h) and spin-flip (-+) configurations. Strong intensity 
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FIG. 5: (Color online) BT7 polarized neutron data measured 
at (0 1.2 0) in the horizontal field configuration, p || Q. (a) 
Comparison of spin-flip (-+) and non-spin-flip (++) scatter- 
ing measured at T = 7 K indicating the magnetic origin of 
the second low energy excitation. Inset: Spin-flip (-+) and 
non-spin-flip (++) scattering of the nuclear reflection (0 2 0). 
(b) Comparison of SPINS high resolution data and the BT7 
subtracted polarized data, the non-spin-flip (++) data was 
subtracted from the spin-flip (-+) data. 



was observed in (++) channel as expected. The observed 
weak (-+) intensity can be attributed to the finite instru- 
mental flipping ratio which we estimate to be ~ 1/9 by 
comparing the integrated intensity of the (-+) and (++) 
scans of the (0 2 0). The spin-flip (-+) and non-spin-flip 
(++) scans at (0 1.2 0) were plotted in Fig. [5] (a) be- 
tween hio — -1.75 meV and 4 meV. Additional magnetic 
scattering was detected in the (-+) spin-flip channel. In 
order to show the peak clearly, the subtracted data, the 
non-spin-flip (++) data was subtracted from the spin- 
flip (-+) data, is plotted together with the SPINS data 
in Fig. 02 (b). At (0 1.2 0), the SPINS data (resolution of 
~ 0.28 meV) shows an excitation centered at 1.37 ± 0.05 
meV. The subtracted BT7 polarized data (resolution of 
~ 1 meV) shows a rather broad peak consistent with the 
SPINS data within experimental error. The polarized 
measurements indicate that this low energy excitation is 
magnetic in origin, which agrees with the temperature de- 
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pendence measurements. However, as shown in Fig. |3l it 
does not fit in the current spin wave model. Attempting 
to analyze the combined dispersions with gaps at ~ 1.2 
meV and ~ 4.7 meV simultaneously using Eq. [2] failed, 
in particular in accounting for the ~ 1.2 meV excitation. 
At this time, the nature of the ~ 1.2 meV excitation is 
not clear, and based on it being nearly dispersionless we 
can only hypothesize that it may be due to a local mag- 
netic excitation. Further studies are necessary in order 
to unravel the nature of this excitation. 



IV. SUMMARY 

The spin dynamics of the ME compound LiCoPC>4 
were determined by inelastic neutron scattering ex- 
periments. Similar to LiNiPO^, LiMnPO^, and 
LiFePC>4 28 , the overall observed magnetic excitation 
spectra in LiCoPCU can be adequately described by a lin- 
ear spin wave theory. Our results indicate that single-ion 
anisotropy is as important as the strong nearest-neighbor 
magnetic coupling and plays an essential role in under- 
standing the spin structure and dynamics of LiCoPC>4. 



However, the observation of the second low energy dis- 
persionless ~ 1.2 meV magnetic excitation is unusual and 
is not contained in the spin wave Hamiltonian, suggest- 
ing that it may be closely related to the strong ME effect 
in LiCoPC>4. The nature of this excitation is not under- 
stood yet and requires further detailed studies. 
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